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, We consider the case of a coupling in the dark cosmological sector, where a dark energy scalar field 

modifies the gravitational attraction between dark matter particles. We find that the strength of 
' the coupling /3 is constrained using current Cosmic Microwave Background (CMB) data, including 

(SJ \ WMAP7 and SPT, to be less than 0.063 (0.11) at 68% (95%) confidence level. Further, we consider 

the additional effect of the CMB-lensing amplitude, curvature, effective number of relativistic species 
and massive neutrinos and show that the bound from current data on /3 is already strong enough 
to be rather stable with respect to any of these variables. The strongest effect is obtained when 
. we allow for massive neutrinos, in which case the bound becomes slightly weaker, /3 < 0.084(0.14). 

A larger value of the effective number of relativistic degrees of freedom favors larger couplings 
between dark matter and dark energy as well as values of the spectral index closer to 1. Adding 
the present constraints on the Hubble constant, as well as from baryon acoustic oscillations and 
supernovae la, we find [3 < 0.050(0.074). In this case we also find an interesting likelihood peak for 
r \ ' /3 — 0.041 (still compatible with at la). This peak comes mostly from a slight difference between 

the Hubble parameter HST result and the WMAP7+SPT best fit. Finally, we show that forecasts 
of Planck+SPT mock data can pin down the coupling to a precision of better than 1% and detect 
Oh whether the marginal peak we find at small non zero coupling is a real effect. 
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I. INTRODUCTION 



Cosmic Microwave Background (CMB) probes have recently broadened our knowledge of primordial acoustic oscilla- 
tions to small angular scales, extending previous measurements of temperature power spectrum (Wilkinson Microwave 
Anisotropy Probe 7, [l[) up to I ~ 3000 with first compelling evidence of CMB lensing from the South Pole Telescope 
(SPT, [2]) and Atacama Cosmology Telescope (ACT, Q). The impact of small-scale CMB measurements and gravi- 
tational lensing on cosmology is relevant [4j and can be used to constrain cosmological parameters and to address one 
of the major issues of present cosmology, that is to say the nature of dark energy |5T-HCl| . The simplest framework for 
' . ■ dark energy models considers dark energy as a cosmological constant A, contributing to about 74% of the total en- 
ergy density in the universe and providing late time cosmic acceleration, while a Cold Dark Matter represents about 



21% (ACDM model). Though theoretically in good agreement with present observations, a cosmological constant 
is somewhat unpleaseantly affected by the coincidence and fine-tuning problems which seem unavoidable in such a 
framework. Many alternative models have been proposed, though it is fair to say that so far no one completely avoids 
these problems. Some encouraging arguments have been put forward in the framework of dynamical dark energy 
models, where a scalar field (quintessence or cosmon) rolls down a suitable potential (lllll2l| possibly interacting with 
dark matter [l3|, [l4| or gravity [HI, [l6| and therefore modifying the growth of structure. Usually, one of the features 
C$ , of such dynamical dark energy models is to have a non-negligible amount of dark energy at early times. The amount 
of early dark energy (early referring to the time of decoupling) influences CMB peaks in various ways and can be 
strongly constrained when including small scale measurements, as shown for instance in Refs. (171 Il8j. 

In this paper, we consider the case of coupled dark energy models, in which dark matter particles feel an interaction, 
additional to gravity, mediated by the dark energy scalar field. Such an interaction introduces effectively a coupling 
between the evolution of the dark energy scalar field and dark matter particles. In this sense, this class of models 
is both an example in which a non-negligible amount of early dark energy is present as well as a typical scenario 
of modified gravity theories. When seen in the Jordan frame, a coupling between matter and dark energy can be 
reformulated in terms of scalar-tensor theories (or f(R) models). This is exactly true when the contribution of baryons 
is neglected. In the Einstein frame, it is common use to neglect a coupling to baryon within coupled dark energy 
models, and consider only dark energy - dark matter interactions. Alternatively, in the Jordan frame, scalar-tensor 
theories (f(R) models) require some sort of screening mechanism (like chameleon (i~9l - [22l | or symmetrons [23]) that 
protects the dark energy scalar field and its mass within high density regions, so that local solar system constraints 
are satisfied. 
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The strength of the coupling affects CMB in several ways, changing the amplitude, the position of the peaks as well 
as contributing to the Late Integrated Sachs- Wolfe (ISW) effect (manifest at large length scales) and to gravitational 
lensing (appearing at small length scales in the temperature spectrum). Moreover, the cou plin g is degenerate with 
the amount of cold dark matter fi c , the spectral index n, the Hubble parameter H(z) (see [24| for a review). After 
recalling the effects of the coupling on CMB, we use a Monte Carlo analysis to constrain the coupling combining 
WMAP and SPT real data. Furthermore we extend our analysis to forecasting the constraints that Planck data are 
expected to put on the coupling parameter, combined with mock SPT data. 

This paper is organized as follows. In Section II we recall the main features of coupled dark energy (CDE) 
cosmologies. In section III we recall effects of the coupling on the CMB spectrum and describe the methods used, 
both with regard to the implementation of the numerical code and the data used for this paper. In Section IV we 
derive the constraints from existing data for several different runs, including effects of the effective relativistic degree 
of freedom N e ff, CMB-lensing, curvature and massive neutrinos. Here we also forecast the constraining capability in 
presence of the forthcoming Planck data, joined with SPT mock data. Finally, in Section V we derive our conclusions. 



II. COUPLED DARK ENERGY 



We consider the case in which an interaction is present between dark energy and dark matter, as illustrated in 
[13, 14, 25-28]. Such cosmologies have to be seen within the framework of modified gravity, since effectively the 
gravitational interaction perceived by dark matter particles is modified with respect to standard General Relativity. 
Coupled dark energy cosmologies considered here are described by the lagrangian: 



C = --d"(pd^ - U{<f>) - m(0)^> + Ad„M , (1) 



in which the mass of matter fields tp is a function of the scalar field </>. We consider the case in which the DE is 
only coupled to CDM (hereafter denoted with a subscript c, while the subscript b will denote baryons). In this case, 
the coupling is not affected by tests on the equivalence principle and solar system constraints and can therefore be 
stronger than the one with baryons. The choice m(<f) specifies the coupling and as a consequence the quantity Q^)^ 
via the expression: 

dlnm(</>) 

= 7T7 PcO^. (2) 



QuAfi ac ts as a source term in the Bianchi identities: 

T u-.n = Q{4>)n ( 3 ) 

If no other species is involved in the coupling then Q( c )/j — ~Q(4>)^- Various choices of couplings have been 
investigated in literature, including constant j3 (29l - [38| and varying couplings [39| . For a constant coupling, typical 
values of fj presently allowed by observations (within current CMB data) are within the range < /3 < 0.06 (at 95% 
CL for a constant coupling and an exponential potential) [2614281 140] . or possibly more [4ll. |42| if neutrinos are taken 
into account or for more realistic time-dependent choices of the coupling. Analysis of the models and constraint on 
these couplings have been obtained in several other ways, including spherical collapse ( [43ll44l and references therein), 
higher-order expansions with the time renormalizazion grou p [45l. TV-body simulations |39l . [jfil It^ and effects on 
supernovae, CMB and cross-correlation of CMB and LSS |27l I28L I40I442L [48l45lT | together with Fisher matrix forecasts 
analysis combining power spectrum and Baryonic Acoustic Oscillations measurements as expected by the Euclid 
satellite (7f| and CMB as expected from Planck (24^ . 

The zero-component of equation ([3]) gives the conservation equations for the energy densities of each species: 

p't = -3Hp<f,(l + w^) - Q WQ , (4) 
p' c = -3Hp c - 



Here we have treated each component as a fluid with T"^ a ^ = (p a + p a )u^u u + p a S^ where = (— a, 0,0, 0) is the 
fluid 4-velocity and w a = Pa/ Pa is the equation of state. Primes denote derivative with respect to conformal time r. 
The class of models considered here corresponds to the choice: 

m(4>) = m e~ l3 ^ T , (5) 
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with the coupling term equal to 



Q(<f>)o = ~J^P c( t > ' ■ ( 6 ) 



Equivalently, the scalar field evolves according to the Klein-Gordon equation, which now includes an extra term that 
depends on CDM energy density: 

0" + 2H0' + a 2 ^ = a 2 /3p c . (7) 
a<p 

Throughout this paper we choose an inverse power law potential defined as: 

v = v^- a (8) 

with a and Vq constants. 



III. THE COUPLING EFFECT ON THE CMB POWER SPECTRUM 

As discussed in (24j, the coupling has two main effects on the CMB: 1) it shifts the position of the acoustic peaks 
to larger £'s due to the increase in the last scattering surface distance (sometimes called projection effect, [35j and 
references therein); 2) it reduces the ratio of baryons to dark matter at decoupling with respect to its present value, 
since coupled dark matter dilute faster than in an uncoupled model. Both effects are clearly visible in Fig. JT]) for 
various values of /3. 

In Figf2]tlie effect of the coupling on the CMB is more evident. The figure shows the quantity Ci{l + l)l 3 /(2ir) (i.e. 
the usual TT spectrum plotted in Fig. ([1) multiplied by I 2 ), as suggested for example in [53]. In Figf2]we also show 
the effect of including CMB-lensing on the TT spectrum for the same value of the coupling constant: the unlensed 
TT spectrum (blue dot dashed) clearly differs from the corresponding lensed one (solid blue) , an effect which is larger 
at small scales (large 1). The position of the peaks remains invariant but the amplitude is larger and, for I 2000, 
the throats appear more pronounced than the peaks. 

The purpose of our analysis is to use recent CMB observations reaching multipoles up to I ~ 3000 to constrain 
dynamical dark energy models in which a coupling is present between the quintessence scalar field (or cosmon, seen 
as the mediator of a fifth force between dark matter particles) and dark matter. In order to do so, we proceed as 
follows. 



A. Theoretical spectra 

Theoretical CMB and lensin g p ower spectra have been produced using the code IDEA (Interacting Dark Energy 
Anisotropics) based on CAMB [54j and able to include dynamical dark energy, Early Dark Energy parameterizations 
(not included in this analysis) as well as interacting dark energy models. In order to include the coupling, both 
background and linear perturbations have been modified following Refs. (30l [35j . The output has been compared to 
an independent code [28| that is built on CMBFAST and the agreement was better than 1%. The difficulty in the 
implementation relies on the fact that the initial conditions cannot be obtained analytically as in simple dark energy 
parameterizations (early dark energy or (u>o, w a )): instead, they must be found by trial and error, through an iterative 
routine that finds the initial conditions required to get the desired present values of the cosmological parameters. 

We have then performed a Monte Carlo analysis integrating IDEA within COSMOMC [55| comparing our theoretical 
predictions with the data presented in the next subsection. 

We recall that the CMB coming from the last scattering surface (LSS) is bent by gravitational structures on the path 
towards us; this effect is called CMB-lensing [J, |56|. The standard deviation of the deflection angle is of the order of 2 
arcminutes, which would correspond to small scales and I > 3000 multipoles, where CMB peaks are already damped 
by photon diffusion. However, deflection angles are correlated with each other over degree scales, so that lensing can 
have an important effect on the scales of the primary acoustic peaks, mainly smoothing them and transferring power 
to larger multipoles. Recently, CMB lensing detection has been claimed by several groups analyzing total intensity 
CMB anisotropies @,|3]. 

CMB lensing is a probe of the expansion rate and naturally depends on the growth of perturbation and on the 
gravitational potentials; since dark energy affects both aspects, CMB lensing represents a way to discriminate among 
dynamical Dark Energy models and ACDM, with promising results [5j4l0(. In particular, the CMB weak lensing 
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Figure 1: CMB TT temperature spectra for three values of /?. Data are taken from WMAP7 [52 



theory in generalized cosmologies has been outlined in [7|. The difference between the lensed and unlensed curves in 
FigfT] shows the typical effects from CMB lensing. The acoustic peaks are smeared because of the correlation between 
different scales induced by lensing, and for the same reason a fraction of power is transferred to the angular domain 
corresponding to the damping tail, therefore dominating that part of the spectrum. Earlier works @ have pointed out 
how the lensing is most relevant in particular in early dark energy models, as it injects power at the onset of cosmic 
acceleration, z ~ 1 ± 0.5 constraining the dark energy abundance in the corresponding epoch. On these lines, the 
analysis in [1C} has shown that the inclusion of lensing data promote the CMB alone to be a probe of the existence 
of dark energy, breaking geometrical degeneracies associated to the pure CMB anisotropics at last scattering. 

Also, the lensing depends on time, combining information from decou plin g (from the last scattering surface of the 
CMB) and z < 5 (when large scale structures formed); recent studies [57ll58fl focus on implementing and investigating 
simulations of CMB lensing through cosmological structures in N-body simulations. 

During matter dominated era (MDE), the potentials encountered along the way are constant in the linear regime 
and the gradient of the potential causes a total deflection angle given by: 



a = -2 



X 'dx fK l x * f , x) v^(xn;T -x), 

Jk(X*) 



(9) 



where x* is the conformal distance of the source acting as a lens, \& is its gravitational potential, r/o ~x is the conformal 
time at which the CMB photon was at position 

One can also define the gravitational lensing potential 



i>{n) ee -2 d X fKiX * *\ *( X n; % - x) 
Jo ]k{x*)jk(x) 



(10) 
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Figure 2: CMB TT temperature spectra (top panel) and dimensionless lensing potential 1(1 + l)C^/(2n) versus multipole I 
for three values of the coupling f3 . 



The lensed CMB temperature Tn in a direction n is given by the unlensed temperature in the deflected direction 
T(h) = T(n') — T(h + a) where at lowest order the deflection angle a — Vi/j is just the gradient of the lensing 
potential. Expanding the lensing potential into spherical harmonics, one can define also the angular power spectrum 
Cf corresponding to the lensing potential, defined as < ipimip*i> m > >— SwS mm tCf; the latter (multiplied by |7(Z + 1)] 2 ) 
is shown in Figj2] (lower panel) for different values of the coupling. As we can see from the plot, the CMB lensing 
potential mainly gives contribution to large scales up to I ~ 1000 or less. However, the lensed CMB temperature 
power spectrum depends on the convolution between the lensing potential and the unlensed temperature spectrum 
(see [4j for more details) whose effect is of several percent at I > 1000, thus being important when estimating the 
spectrum up to small scales of / ~ 3000, as for the data we consider in the following. 
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Figure 3: WMAP7 and SPT data used for this work. SZ effects are not included in the best fit line (solid black) (see fig. 5 of 
3 for the best fit including SZ). The best fit from WMAP+SPT, for a ACDM without foregrounds is also shown (dash-dotted 
black line). This plot is similar to fig. 5 of Q (without SZ in the best fit); we reproduced it here for convenience and because 
it will allow us to neglect SZ when combining Planck with mock SPT data, in section V. All foregrounds are instead included 
when SPT is combined with WMAP7. 



B. Observed data 



We have compared theoretical predictions of CMB lensed spectra with two datasets. The first includes WMAP7 
temperature spectral]]] . The second includes the recently released power spectrum data from SPT [2]. Together with 
the ATC [H, SPT [2|| has recently shown evidence of CMB-lensing with an enhancement in the CMB temperature 
power spectrum up to I ~ 3000. We do not use ACT data in this analysis. 

In order to use SPT data, we have installed the likelihood provided by the SPT team [2] on their SPT website 
[7rj| and integrated it with the recommended version of COSMOMC [77} (August 2011). We have then implemented 
IDEA on this version. Care has to be used whenever small multipoles in the range 2000 — 3000 are used, due to 
several sources of foregrounds active on those scales. In particular, whenever SPT data are included in the analysis 
we also marginalize over the three nuisance parameters described in Q: two of them refer to Poisson point sources 
and clustered point sources; the third one adds power from the thermal and kinetic Sunyaev-Zel'dovich (SZ) effects, 
which are also an example of secondary CMB anisotropics. These effects are relevant when small scales (I £ 2000) 

are included and have to be taken into account whenever SPT data are used. In Fig.© and Fig. J2J) we show the 
superposition of WMAP7 and SPT data, together with the best fit of the combined set. 

The baseline set of parameters includes = Vl^h 2 ,il c h 2 ,fld e ,0s,logA,n s ,T. As in other papers, these parameters 
represent the fractional abundances of the various species, as well as the amplitude and shape of the primordial power 
spectrum, and the reionization optical depth; as stressed already, we also consider the three nuisance parameters when 
SPT is included; when we impose spatial flatness, the present dark energy density fide becomes a derived parameter; 
in addition, when CQ is included two more parameters are added: (3 and a. Again, j3 represents the coupling between 
dark matter particles in eqfSl a is the parameter in the scalar field potential ((5J) that drives the long range interaction. 
Different runs are illustrated in Tab|H The Helium abundance Yjie is derived following BBN consistency (see for 
details). 
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Figure 4: WMAP7 and SPT data used for this work, with an extra factor of £ 2 . SZ effects are not included in the best fit line 
(solid black) (see fig.5 of for the best fit including SZ). The best fit from WMAP+SPT, for a ACDM without foregrounds 
is also shown (dash-dotted black line). 



Run 


CMB-Lensing 


WMAP7 


SPT 


Planck 


Parameters 


cql 


/ 


/ 


/ 


X 


baseline + a + f3 


cq2 






X 


X 


baseline + a + f3 


cq4 






/ 


X 


baseline + a + f3 + 7V„ 


cqlNL 


X 






X 


baseline + a + f3 


cqlK 


/ 


/ 


/ 


X 


baseline + a + j3 + curvature 


cqlAL 


/ 






X 


baseline + a + j3 + Ai 


cqlv 


/ 






X 


baseline + a + f3 + f v 


cqlhst 


/ 






X 


baseline + a + ft + HST + BAO + SNae 


cql PI 


/ 


X 


/ (mock) 


/ (mock) 


baseline + a + /3 + HST + BAO + SNae 



Table 1: COSMOMC Monte Carlo simulation runs described in this paper; they all refer to CQ models. 



IV. RESULTS 

As a first step, we have performed a run using WMAP7 and SPT with a ACDM model, and we find results 
compatible with Q. Note for the following, that in [2j the authors report the mean values of each parameter, together 
with its standard deviation. We instead report the best fit values and the marginalized errors at 68% and 95% 
confidence level. We now describe results from the runs illustrated in Tab|H 

A. Baseline plus /3 and a 

The first two Monte Carlo runs we describe (cql and cq2) use the baseline set of parameters O = 
{flbh 2 , fl c h 2 , 9 S , logA, n Sl r}; in addition, two more parameters are added to account for the coupling (/?) and the 



ACDM best fit + PS + CL 
ACDM best fit 




Best fit values for coupled quintessence 



Parameter 


cql (WMAP7 + SPT) 


cq2 (WMAP7) 


Q b h 2 


n n99 +0.0007 
UiU "-0 .00013 


n n9s+ 00044 

u - uzo -0. 00070 


tt c h 2 


n i i +0.0022 
U.ll_o.oil 


n -i -1+0.0019 
U - XJ --0.016 


6 S 


, n4 +0.0024 
L - u ^-0. 00072 


, n4 +0.0027 
1 - u ^0.0025 


r 


u.u»i_ 012 


0-089 .0072 


n s 


q fi +o.oi9 

u - 3u -0.0056 


07+0.021 

u - y ' -0.012 


w 


-0.88™ 


-0 Q7+ 017 
u - y '-0.03 


P 


n ni2 +0 050 

u.uiz_ 012 


0.0066?:8Sb86 


P 


< 0.063(0.11) 


< 0.078(0.14) 


a 


ri 99+O.28 
u - zz -0.090 


n -,0+0.37 

u - lo -0.004g 


fide 


n 79 +0.076 
u - ' z -0.012 


n 7 9+0.093 
u - ' z -0.016 


Age/Gyr 


1 ■? Q+0.012 

10.0_g 39 


13.818:? 


z re 


10.9±?;! 3 




Ho 


68.7to'gg 


69.3li° 7 7 


n sz 

-^3000 






-^3000 






n CL 
^3000 


r Q + 1.9 
°-°—2.4 




-Log (Like) 


3756 


3737 



Table II: Best fit values and 1-er errors comparing runs cql and cq2. Both runs include coupling; cql uses WMAP7 + SPT 
while cq2 uses WMAP7 only. For f) we also write in brackets the value of the 2a marginalized error. 

dark energy scalar field potential (a). Run cql compares theoretical spectra with WMAP7+SPT data; run cq2 in- 
cludes WMAP7 data only. Results are shown in TabOH where we report best fit values with 1-sigma (68%) errors on 
various parameters. For /3 we also report upper 1 and 2 a limits. When WMAP7 only is considered, the coupling is 
constrained to be /3 < 0.078(0.14) at 1 (2) a, while including SPT data, and therefore small scales and large multipoles, 
the coupling is constrained down to < 0.063(0.11). 

In the past, [13 found fj A < 0.16 at 95% c.l. (note that our definition of /3 is (3 = yf{2/$)P A = 0.13 at 2 a) 
using COBE, Boomerang, Maxima, DASI and fixing the optical depth r. [28[ found j3 < 0.061(0.11) at 1 (2) a using 
WMAP1 data and < 0.11(0.16) for pre-WMAP data (using a different set than (27j). Ref. found C < 0.034 
(0.066) (in their notation C = (3 A = y/V^P so that their result is equivalent to (3 < 0.028(0.054)) using WMAP5, 
SNLS, HST, LRG, SDSS. See also [HO] (with a different definition of the coupling) and [4l| (where massive neutrinos 
were included). These constraints are not easily comparable to ours since in these papers different priors have been 
used and/or some parameters have been kept fixed. Overall, however, we find an agreement with the more recent 
constraints to within a factor of 50% at most. 

The 2D confidence contours are plotted in FigJS] Here we show a selection of the most interesting likelihood 
contours vs the coupling 0. In FigJ5]we also show ID likelihood contours, comparing results from WMAP7+ SPT 
with WMAP7 only. Note that there is no dependence of cosmological parameters from a, as expected since a only 
affects late time cosmology; the range in a was therefore safely chosen to be between 0.13 and 0.5, small enough to 
get reasonable speed for the runs. The value of w is arbitrary and approximately related to a via the expression: 
w = — 2/ (a + 2); the interval chosen for a is such that w still assumes reasonable values, at least smaller than -0.8. 
As it appears clearly from Fig. [HI CMB is practically insensitive to a or w within the range we consider. 



B. Effective number of relativistic species N e ff 

The effective number of relativistic species before recombination, usually denoted by N e ff is higher than the 
number of relativistic neutrino species (3.046) due to photons produced in electron-positron annihilation at the end 
of neutrino freeze-out [H-[62|. Using_WMAP7 + SPT data for a ACDM, N eff = 3.046 was found to be preferred over 
zero relativistic species (N e ff =0) (2J. If N e ff is left free to vary and marginalized over, its best fit value can be even 
larger: found N eff > 2.7 at 95% CL using WMAP7 alone; (63j found N eff = 5.3 ± 1.3 using WMAP7 + ACT; Q 
found N e ff — 3.85 ±0.62 using WMAP7 + SPT. If relativistic species are present, the expansion rate during radiation 
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Figure 5: Confidence contours for the cosmological parameters for coupled quintessence models. We compare runs cql (red) 
and cq2 (white). The light blue asterisks mark the best fit points for cql. 1-sigma and 2-sigma contours are shown. 



dominated era increases |2|. I64j - l66j | . We have redone the analysis in the case of coupled quintessence, to check whether 
the number of relativistic degrees of freedom is degenerate with the coupling. The effect of marginalizing over N e ff 
on the coupling is shown in Tab lHIl In FigJT] we plot the likelihood contours for a selection of parameters, vs the 
coupling /3, comparing different runs. 

Contours are larger when we allow N e ff to vary, but the range in f3 is not affected considerably ((3 < 0.074(0.12) 
instead of /3 < 0.063(0.11)). When a coupling is present, we find that the best fit for the number of relativistic species 
is given by N e ff — 3.84^q'4q when using WMAP and SPT, similar to the value mentioned before and evaluated in 
absence of a coupling. The allowed range for Yu e , obtained through BBN consistency, increases a lot when N e ff is 
free to vary (65[. N e ff is degenerate with dark matter and the spectral index, which in turn are weakly degenerate 
with /3, though no direct degeneracy appears between j3 and N e ff, as shown in FigjS] 

It is interesting to see (Fig(7j) that when we allow for an effective number of relativistic degrees of freedom, marginal- 
izing over N e ff, the coupling from WMAP7+SPT data increases to a best fit value of f3 ~ 0.03. A larger value of 
N e ff (best fit ~ 3.8) favors larger couplings between dark matter and dark energy as well as values of the spectral 
index closer to -1 (n s ~ 0.99). 
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Figure 6: ID likelihoods for the cosmological parameters for coupled quintessence models. We compare runs cql (solid blue) 
and cq2 (dashed red). As expected, in coupled quintessence models, no dependence is seen on the value of a and w, which are 
arbitrary. We restrict the analysis to a reasonable range in w and a range in a small enough to speed up the iterative routine 
that finds initial conditions. 



C. CMB lensing 

In order to test the effect of CMB-lensing on coupled quintessence, we have redone a run cql (WMAP + SPT) 
without lensing. The presence of a constant coupling doesn't seem to be very much affected by lensing in the TT 
CMB spectra, as we can see in Fig.® where we compare run cql with run cqlNL. If no lensing is included, the bound 
on (3 is slightly (but not significantly) larger: (3 < 0.068(0.13). 
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Best fit values for coupled quintessence WMAP7 + SPT 



Parameter 


cql (baseline + /3+a) 


cq4 (baseline + (3+a+N e ff) 


n b h 2 


n n99 +0.0007 
UiU "-0 .00013 


n n9^+ - 00054 
u.uzo_ 000 54 


tl c h 2 


n i i +0.0022 
U.ll_o.oil 


n 19+0.0076 
u - lz -0.016 


Os 


, n4 +0.0024 
L - u ^-0. 00072 


i n4 +0.0012 
i - uq; -0.0025 


T 


u.u»i_ 012 


0.086^0.0^30 


n s 


n OR+0 019 

u.»o_ 0056 


n qq+0 029 
u.ya_ .oi4 


w 
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Table III: Best fit values and 1-a errors comparing runs cql and cq4- Both runs include coupling and use WMAP7 + SPT; in 
addition, cq4 marginalizes over N e ff. For f3 we also write in brackets the value of the 2a marginalized error. 

Furthermore, similarly to one can rescale the lensing potential power spectrum by a factor At,: 

<3f* -> A L Cf* (11) 

All runs discussed so far fix At, = 1. In order to test the effect of lensing we also performed a run (cqlAL) in which 
we vary Al and marginalize over it. The Aj parameter was found to be At, = 0.94 ± 0.15 when using WMAP + 
SPT with a ACDM model [lH (see also 0,(67- 69] for different datasets). When a coupling between dark matter and 
dark energy is included, we find that the best fit for Al is At, = 0.861q'?2j still compatible with one. The bound on 
the coupling is of the same order as in the case in which At, is fixed: j3 < 0.063(0.11). In other words, we don't gain 
much marginalizing over At, instead of fixing it to one, given that Al best fit is very close and fully compatible with 
one. Though with the data considered here the effect is not significant, At, is also correlated with dark matter and 
n s , which in turn are correlated with j3, as shown in fig llOl 

D. Curvature 

If we release the constraint of a flat universe and allow for curvature and coupling (run cqlK), we get = 
— 0.0068lQ Q3g 2 , which is compatible with a flat Universe. In this case, the constraint on /3 is slightly less restrictive, 
< 0.071(0.13) but the bound on (3 is already stringent enough not to be affected so much by the uncertainty on 
curvature. Contours are shown in Figll II where they are compared to run cql, in which a flat universe was assumed. 
We also show in Fig[l2] how curvature is degenerate, as expected, with the Hubble parameter, though no direct 
degeneracy is seen between SIk and the coupling /3. 

E. Massive neutrinos 

Up to now, we have fixed the relative fraction of massive neutrinos /„ to zero. We now consider run cql v in which 
we also allow for a non zero fraction of massive neutrinos and marginalize over /„. In this case the range allowed for 
the coupling is /3 < 0.084(0.14), slightly bigger than when using massless neutrinos, as expected |4]J. We update the 



12 



1.05 




Figure 7: Confidence contours for the cosmological parameters for coupled quintessence models. We compare runs cql (blue), 
cq2 (green) and cq4 (yellow). The light blue asterisks mark the best fit points for cql while the pink asterisks mark the best 
fit points for cq4. 1-sigma and 2-sigma contours are shown. 

results of [41] using both WMAP7 and SPT data. The degeneracy between massive neutrinos and /3 is clearly shown 
in FigQJ The best fit value for /„ is f v = 0.065±£^ 



F. Combining WMAP7 and SPT with HST, BAO and Supernovae la data 

As discussed earlier on in this paper and in [24j|, the coupling /3 is degenerate with the Hubble parameter. In order 
to investigate the effect of this degeneracy on the constraints from data, we did another run (cqlhst) in which we 
combined the data used for cql (WMAP7 and SPT) to also included baryon acoustic oscillations (BAO) (7^1, Hubble 
Space Telescope constraints on H (HST) [7l| and Supernovae la (SNae) data [zl] as from COSMOMC (Aug 2011). 
Likelihood contours for this case are shown in FigO Best fit values and errors are shown in Tab II VI left column. 
There seems to be an interesting preference for a non zero coupling, though the values are clearly still compatible with 
zero at la. This peak comes mostly from a slight tension between the Hubble parameter HST result (ft = 0.738±0.024) 
and our WMAP7+SPT best fit for = (ft = 0.685 ± 0.025 ). Notice however that even for /3 = we are not in 
an exact ACDM since in our model w is close, but not exactly equal, to -1. It is interesting then to test whether the 
forecoming data from Planck can confirm or reject this non zero coupling. This we do in the next section. 
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N eff 

Figure 8: Confidence contours for a choice of cosmological parameters vs N e f / within coupled quintessence models for run cq4- 
1-sigma and 2-sigma contours are shown.. AT e // is degenerate with floAih 2 , n a , which are in turn degenerate with /3. No direct 
degeneracy can however be seen between j3 and N e ff. The pink asterisks mark the best fit for cq4. 

G. Combining Planck and SPT mock data 

As a further analysis, we forecast the effect that Planck data would have on the coupling parameter, when combined 
with the power measurement of SPT (run cqlPl). Since Planck data are not yet available, we produce a set of mock 
data [73]. We have therefore implemented FutureCMB [74[ in our modified version of COSMOMC and, using as 
fiducial power spectrum a ACDM model Df 1 , with D\ — 1(1 + 1)/2ttCi (black dot-dashed in FigfT5|. We have then 
produced an SPT mock spectrum as follows: 
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Figure 9: Likelihood contours for cosmological parameters in presence of a coupling for WMAP+SPT (run cql, including 
lensing, red contours), as compared to the same run done without CMB lensing (white contours). Light blue asterisks mark 
the best fit points for cql. 

we have added to the same fiducial model used to generate Planck mock data, the effect of poisson sources (PS) 
(blue dotted) using the value for the nuisance parameters Di PS 3000 = l8.1^iK 2 with a dependence from the 
momentum oc (Z/3000) 2 @|; 

we have added to D\ h + PS the effect of clustered sources (CL) (light blue dashed line), using the value for the 
nuisance parameters Di CL 3Q0Q — 3.5fj,K 2 with a dependence from the momentum oc (Z/3000)°.8 for all I. This is 
not entirely correct as the dependence is slightly different for I > 1500 and I < 1500 but differences are thought 
to be small @; we obtain the D t = Dj h + PS + CL. We neglect here the effect of SZ. 

we have then convolved the L>i for the SPT window functions and created an SPT mock data with Di given by the 
convolved Di and errors given by SPT data from The final SPT mock data is plotted in Fig fT5|) (red dotted 
line) . 



For this case, we also include BAO, HST and Supernovae la data as from Cosmomc (August 2011 version), to 
break the degeneracy between the coupling and the Hubble parameter. Some representative 2D confidence regions 
are in Fig. [16] We then find that (3 < 0.012(0.030); though stronger than WMAP7+SPT, this is still a pessimistic 
bound, since it includes Planck but still considers SPT data with errors released by [H; by the time Planck data will 
be available, better SPT (or ACT) data may have been released. In Table HVl we report the best fit values, together 
with the left and right errors at 68% and 95% CL around the best fit. The left column refers to run cqlhst, done with 
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Figure 10: Likelihood contours for cosmological parameters when Al is allowed to vary. 

real data (WMAP+SPT+HST+BAO+SNae); the right column shows the forecasted values for Planck+SPT mock 
data, plus SNae, HST and BAO. 

We finally plot in Fig. [XT] a comparison between current data and future observations, marginalizing over all 
parameters except /3. It is clear that Planck data will reach a precision sufficient to tell whether the peak in ft is a 
real detection or just a fluke. 

V. CONCLUSIONS 

We have considered the possibility that the evolution of dark matter and dark energy might be connected by a 
constant coupling, of the type illustrated in (l3l . We have used current CMB data from WMAP7 and SPT to 
constrain the coupling parameter j3. We find that j3 is constrained to be less than 0.063 (0.11) at 68% (95%) C.L. 
when SPT data are included, with respect to j3 < 0.078(0.14) coming from WMAP7 only. We have done a number 
of tests to check whether this bound depends on the degeneracy with other parameters (lensing, curvature, massive 
neutrinos, N e ff, HST/BAO/SNae data). If the effective number of relativistic degrees of freedom N e ff is allowed 
to vary, no much gain is obtained on f3, which still needs to be j3 < 0.074(0.12). We have further considered the 
effect of CMB-lensing, both with a run which includes no lensing and by marginalizing over Al, a parameter which 
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Figure 11: Likelihood contours for cosmological parameters of run cql (red contours) as compared to run cql, in which SIk is 
allowed to vary (white contours). Light blue asterisks mark the best fit points of cql while yellow asterisks mark the best fit 
points for cqlK. 



encodes the rescaling of the lensing power spectrum. Al is slightly degenerate with n s ,QDMh 2 which in turn are 
degenerate with /3, though no direct degeneracy is seen between f3 and Al- If the assumption of a flat universe is 
released, constraints on j3 weaken back almost to the level of constraints given by WMAP only (flat universe), with 
/3 < 0.071(0.13). Degeneracy with massive neutrinos widens the coupling constraints to be (3 < 0.084(0.14) when we 
marginalize over the fraction of massive neutrino species /„. We conclude that the bound on j3 from current data 
is already strong enough to be quite stable with respect to a better knowledge of other parameters and to all cases 
considered. 

When WMAP+SPT are considered (run cql), the best fit value for /?, though still fully compatible with zero, has 
a best fit of j3 = 0.012+q "12- It i s interesting to see that when we allow for an effective number of relativistic degrees 
of freedom, marginalizing over N e ff the coupling from WMAP7+SPT data increases to a best fit value of j3 ~ 0.03. 
A larger value of N e f / favors larger couplings between dark matter and dark energy and values of the spectral index 
closer to 1. Including SPT data does not improve significantly constraints on the coupling f3. Inclusion of additional 
priors from HST, BAO and SNae moves the best fit to f3 = 0.041, again still compatible with zero at la. We forecast 
that the inclusion of Planck data will be able to pin down the coupling to about 1% and therefore detect whether the 
small non-zero coupling present in current data is washed away with more data. 
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Figure 12: Likelihood contours for cosmological parameters when Q,k is allowed to vary. Yellow asterisks mark best fit points 
for run cqlK. 



Acknowledgments 

V.P. is supported by the Marie Curie IEF, Project DEMO - Dark Energy Models and Observations. Support was 
given to V.P. and C.B. by the Italian Space Agency through the ASI contracts Euclid-IC (1/031/10/0). V.P. thanks 
Joanna Dunkley for useful discussion and advices. We thank Christof Wetterich for helpful comments. C.B. also 
acknowledges support from the PD51 INFN initiative. LA. acknowledges support from DFG - Deutsche Forschungs- 
gemeinschaft through the project TRR33 "The Dark Universe". 



[1] E. Komatsu, K. M. Smith, J. Dunkley, C. L. Bennett, B. Gold, G. Hinsh aw, N. Jarosik, D. Larson, M. R . Nolta, L. Page, 

et al, The Astrophysical Journal Supplement Series 192, 57 (2010) , [arXiv : 1001 .T5 38v3 [astro-ph . C0]~1 

[2] R. Keisler, C. Reichardt, K. Aird, B. Benson, L. Bleem, et al., Astrophvs.J. 743. 28 (20111, larXiv: 1105 . 31821 

[3] S. Das, B. D. Sherwin, P. Aguirre, J. W. Appel, J. Bond, e t al, Phys.Rev.Lett. 107, 021301 (2011), larXiv: 1103.21241 

[4] A. Lewis and A. Challinor, Physics Reports 429, 1 (2006) , [arXiv : arXiv : astro-ph/0601594 1 



18 



0.14 




c 1 " 0.96 




0.05 0.1 0.15 0.2 
f 

v 



Figure 13: Likelihood contours for cosmological parameters when f u is allowed to vary. Pink asterisks mark best fit points for 
run cq\v. 



L. Verde and D. N. Spergel, Phys. Rev. D 65, 043007 (2002), URL http : //link . aps . org/doi/10 . 1103/PhysRevD . 65 . 043007 

F. Giovi, C. Baccigalupi, and F. Perrotta, ArXiv Astrophysics e-prints (2003), arXiv: arXiv: astro-ph/0309422 

V. Acquaviva, C. Baccigalupi, and F. Perrotta, Phys. Rev. D 70, 023515 (2004), arXiv:arXiv:astro-ph/0403654 

V. Acquaviva and C. Baccigalupi, Phys. Rev. D 74, 103510 (2006), [arXiv : arXiv : astro- ph/05076 44| ~ 

W. Hu, D. Huterer, and K. M. Smith, Astrophys. J. Lett. 650, L13 (2006) , |arXiv : arXiv : astro-pn /0607316 

B. D. Sherwin, J. Dunkley, S. Das, J. W. Appel, J. Bond, et al., Phys.Rev.Lett. 107, 021302 (2011), ar Xiv: 1105.04191 

C. Wetterich, Nucl. Phys. B302, 668 (1988). 

B. Ratra and P. J. E. Peebles, Phys. Rev. D3 7, 3406 (1988). 

L. Amendola, Phys. Rev. D62, 043511 (2000), |arXiv: astro-ph/9908023"] 

V. Pettorino and C. Baccigalupi, Phys. Rev. D77, 103003 (2008). larXiv: 0802 . 10861 

S. Matarrese, C. Baccigalupi, and F. Perrotta, Phys. Rev. D70, 06 1301 (2004), |arXiv: astro-ph/04 03480 
F. Perrotta and C. Baccigalupi, Phys. Rev. D65, 123505 (2002), [arXiv : astro-ph/020 1335 

E. Calabrese, R. de Putter, D. Huterer, E. V. Linder, and A. Melchiorri, Phys. Rev. D 83, 023011 (2011), 
larXiv: 1010.56121 

C. L. Reichardt, R. de Putter, O. Zahn, and Z. Hou (2011). larXiv: 1110.53281 

J. Khoury and A. Weltman, Phys. Rev. D 69, 044026 (2004) , [arXiv : arXiv : astro-ph/0309411 1 
L. Hui, A. Nicolis, and C. Stubbs, Phys. Rev. D80, 104002 (2009), arXiv: 0905. 2966 



19 



1.046 




Figure 14: Likelihood contours for cosmological parameters when data from HST, BAO and Supernovae la are included. Blue 
asterisks mark best fit points for run cqlhst. 



A. Upadhye, W. Hu, and J. Khoury, ArXiv e-prints (2012) . larXiv : 1204 . 39061 

A.-C. Davis, E. A. Lim, J. Sakstein, and D. J. Shaw, ArXiv e-prints (2011) , larXiv : 1 102 . 52781 

K. Hinterbichler, J. Khoury, A. Levy, and A. Matas, Phys. Rev. D 84, 103521 (2011). larXiv: 1107 . 2fT2"l 

L. Amendola, V. Pettorino, C. Quercellini, and A. Vollmer, Phys. Rev. D 85, 103008 (2012), larTiv: lill ■ 14541 

H. Kodama and M. Sasaki, Prog. Theor. Phys. Suppl. 78, 1 (1984). 

L. Amendola, Phys. Rev. D69, 103524 (2004), [arXiv: astro-ph/0311175| 

L. Amendola, C. Quercellini, D. Tocchini-Valentini, and A. Pasqui, Astrophys. J. Lett. 583, L53 (2003), 

arXi v:astro^ ph/0205097| 

L. Amendola and C. Quercellini, Phys. Rev. D 6 8, 0235 14 (2003 ), |arXiv: astro-ph/030 3228 
C. Wetterich, Astron. Astrophys. 301, 321 (1995), |arXiy? hep-th/9408025l 

L. Amendola, Physical Review D 62, 43511 (1999) JarXT7: 9908023 

G. Mangano, G. Miele, and V. Pettorino, Mod. Phys. Lett. A18, 831 (2003), [arXiv: astro-ph/0212518"| 
L. Amendola, Physical Review D 69, 103524 (2004) , larXiv : 031 1 175v2"l 

T. Koivisto, Phys. Rev. D72, 043516 (2005) , [arXiv : astro- ph/0504571 1 

Z.-K. Guo, N. Ohta, and S. Tsujikawa, Phys. Rev. D76, 023508l 2007JjarXiv:astro-ph/0702015| 
V. Pettorino and C. Baccigalupi, Phys. Rev. D77, 103003 ( 2008) . larXiv : 0802 . 10^61 

C. Quercellini, M. Bruni, A. Balbi, and D. Pietrobon, Phys. Rev. D 78, 063527 (2008) . larXiv : 0803 ■ 19761 
M. Quartin, M. O. Calvao, S. E. Joras, R. R. R. Reis, and I. Waga, JCAP 0805, 007 (2008), arXiv: 0802. 0546 



20 




Figure 15: SPT mock data as obtained from the fiducial ACDM spectrum. 



[38] J. Valiviita, R. Maartens, and E. Majerotto, Mon.Not.Roy.Astron.Soc. 402, 2355 (2010), [arXlv; 0907 . 4987 
[39] M. Baldi, MNRAS 411, 1077 (2011), larXiv : 1005 ■ 21881 

[40] R. Bean, E. E. Flanagan, I. Laszlo, and M. Trodden, Phys. Rev. D78, 123514 (2008) . larXiv : 0808 ■ llg5l 
[41] G. La Vacca, J. R. Kristiansen, L. P. L. Colombo, R. Mainini, and S. A. Bonometto, JCAP 0904, 007 (2009), 
larXiv: 0902 .27111 

[42] J. R. Kristiansen, G. La Vacca, L. P. L. Colombo, R. Mainini, and S. A. Bonometto, New Astron. 15, 609 (2010), 
larXiv: 0902 .27371 

[43] N. Wintergerst and V. Pettorino, Phys.Rev. D82, 103516 (2010), ar Xiv: 1005. 12781 

[44] R. Mainini and S. Bonometto, Phys. Rev. D74, 043504 (2006), |arXiv:astro-ph/060562l"1 

[45] F. Saracco, M. Pietroni, N. Tetradis, V. Pettorino, and G. Robbers, Phys.Rev. D82, 023528 (2010) . larXiv : 091 1 . 53961 

[46] M. Baldi and V. Pettorino, Mon.Not.Roy.Astron.Soc. 412, LI (2011). larXiv : 1006 ■ 37BT1 

[47] M. Baldi, V. Pettorino, G. Robbers, and V. Springel, Mon. Not. Roy. Astron. Soc. 403, 1684 (2010). 

[48] R. Mainini and D. F. Mota (2010), arXiv: 1011 . 00831 

[49] F. De Bernardis, M. Martinelli, A. Melchiorri, O. Mena, and A. Cooray (2011), larXiv: 1104. 06521 
[50] J.-Q. Xia, Phys.Rev. D80, 103514 ( 2009) , larXi v : 09 1 1 . 4821)1 

[51] M. Martinelli, L. L. Honorez, A. Melchiorri, and O. Mena, Phys. Rev. D 81, 103534 (2010) , larXiv : 1004 . 24101 

[52] D. Larson, J. Dunkley, G. Hinshaw, E. Komatsu, M. Nolta, et al., Astrophys.J.Suppl. 192, 16 (2011), larXiv : 1001 . 4635 

[53] R. Hlozek, J. Dunkley, G. Addison, J. W. Appel, J. Bond, et al. (2011), arXiv: 1105.4 88T1 

[54] A. Lewis, A. Challinor, and A. Lasenby, Astrophys.J. 538, 473 (2000), arXiv:astro-ph/9911177 

[55] A. Lewis and S. Bridle, Phys. Rev. D 66, 103511 (2002) , [arXiv : arXiv : astro -ph/0205436 

[56] M. Bartelmann and P. Schneider, Phys. Rept. 340, 291 (2001), |arXiv:astro- ph/9912508| 

[57] C. Carbone, C. Baccigalupi, M. Bartelmann, S. Matarrese, and V. Springel, MNRAS 396, 668 (2009) . larXiv : 0810 . 4f4"5l 
[58] R. Teyssier, S. Pires, S. Prunet, D. Aubert, C. Pichon, A. Amara, K. Benabed, S. Colombi, A. Refregier, and J.-L. Starck, 
Astronomy and Astrophysics 497, 335 (2009) , larXiv : 0807 . 36BT1 



21 



Best fit values for coupled quintessence 
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Table IV: Best fit values and 1-a errors comparing runs cqlhst and mock cqlPl including SPT 
SN data. For /3 we also write in brackets the value of the 2a marginalized error. 
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Figure 16: Confidence contours for Planck + SPT mock. Blue asterisks mark best fit points for run cqlPl. 



22 



cq1 

cqlPI 

cqlhst 




Figure 17: ID likelihood for the coupling /3 from runs cql, cqlhst, cqlPI. We recall that cgl and cqlhst are based on real data 
while cqlPI estimates the forecasted constraints from mock data, around a fiducial ACDM cosmology. 
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